Varying couplings in the early universe: correlated variations of a and G 
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The Cosmic Microwave Background anisotropics provide a unique opportunity to constrain simul- 
taneous variations of the fine-structure constant a and Newton's gravitational constant G. Those 
correlated variations are possible in a wide class of theoretical models. In this brief paper we show 
that the current data, assuming that particle masses are constant, gives no clear indication for such 
variations, but already prefers that any relative variations in a should be of the same sign of those 
of G for variations of ~ 1%. We also show that a cosmic complementarity is present with Big 
Bang Nucleosynthesis and that a combination of current CMB and BBN data strongly constraints 
simultaneous variations in ct and G. We finally discuss the future bounds achievable by the Planck 
satellite mission. 
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I. INTRODUCTION 

The behavior of nature's fundamental couplings is a 
subject of much recent interest. There is ample ex- 
perimental evidence showing that they run with energy, 
and many particle physics and cosmology models sug- 
gest that they should also roll with time. This explains 
why the European Space Agency (ESA) and the Euro- 
pean Southern Observatory (ESO) now list varying fun- 
damental constants among their key science drivers for 
the next generation of facilities. Recent technological de- 
velopments have for the first time provided us with tools 
to accurately test this hypothesis. Two recent reviews 
Jl], H discuss most of these developments. 

In particular, there are controversial claims for time 
variations of the fine-structure constant a || and the 
proton-to-clcctron mass ratio fi j|| at rcdshifts z ~ 1 — 3, 
as well as of possible spatial variations of the latter in the 
Galaxy || . While these have so far not been confirmed by 
independent analyses [^-^j], their potential implications 
certainly warrant further study. In the present work we 
focus on the cosmic microwave background (CMB) as a 
means to probe the behavior of these couplings at red- 
sifts z ~ 10 3 , but we will also discuss the importance of 
constraints at much higher rcdshifts (z ~ 10 10 ), coming 
from Big Bang Nucleosynthesis. 

Typically, in any sensible theory where a coupling is 
rolling, one generically expects the others to do so as well, 
though possibly at fairly different rates. Since the rolling 
is expected to be due to the same underlying mechanism 
(the most natural of which will be a dynamical, funda- 
mental scalar field), the rates of change of the various 
couplings will in fact be related in any given theory 

From an experimental point of view, one can either 
take the simplifying assumption that only the coupling 



one wants to constrain is varying while the others are con- 
stant, or try to constrain joint variations at the expense of 
a more complicated analysis. In the latter case one can ei- 
ther choose a particular set of such relations (thereby con- 
straining only a particular theory, though usually quite 
tightly) or phenomenologically treat the different varia- 
tions as independent (thus obtaining model-independent 
but usually weaker constraints). In this work we will 
use the CMB to constrain possible variations of a and 
Newton's constant G, assuming them to be related by a 
phenomenological parameter whose value will be differ- 
ent in various fundamental physics scenarios. 

The importance of a as a fundamental physics probe 
stems from the fact that it is ubiquitous in electromag- 
netic processes. In the past this has been extensively used 
to constrain the fine-structure constant [jl3] pOf . These 
studies yielded results consistent with no variation, but 
due to degeneracies with other cosmological parameters 
the accuracy is lower than that of low-redshift measure- 
ments: it is only with the latest available data that con- 
straints stronger than the percent level have been ob- 
tained 

Speaking of variations of dimensional constants obvi- 
ously has no physical significance: one can design any 
variation by defining appropriate units of length, time 
and energy. However, one is free to choose an arbitrary 
dimensionful unit as a standard and compare it with 
other quantities. If one assumes particle masses to be 
constant (as we will do in this paper), constraints on the 
gravitational constant G are in fact constraining the di- 
mensionless product of G and the nucleon mass squared. 
With this caveat, constraints on a rolling G provide key 
information on the gravitational sector. Somewhat para- 
doxically, G was the first constant to be measured but 
is now the least well known, a result of the weakness of 
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gravity. Indeed, in the past two decades our knowledge of 
its value didn't substantially improve from the precision 
of 0.05% reached in 1942 (sec |22 ). Recent laboratory 
measurements (see e.g. [|23|) point towards an uncer- 
tainty at the level of ~ 0.4%, while other works claim an 
improved precisions below 0.01% (|24|). Analysis of the 
secular variation of the period of nonradial pulsations of 
the white dwarf G117-B15A (|2[|) has produced comple- 
mentary constraints at ~ 0.1% level. 

On the other hand, as shown in p(| a variation in 
the gravitational constant could also affect the CMB 
anisotropy spectra, with current bounds of the order of 
~ 10% (sc e.g. j27J) and constraints at level of ~ 1% 
achievable with future CMB experiments as Planck. 

The interesting point is that the CMB is an observ- 
able potentially sensitive to variations in both fundamen- 
tal constants. It is therefore timely to perform a com- 
bined analysis of CMB data considering simultaneous 
variations in a and G in order to investigate the possible 
correlations and deviations from the standard values. 

Specifically, we will consider that the variations of a 
and G are related by 



Aa _ AG 
a G 



(1) 



with Q a free parameter that can be positive or negative, 
but not much larger than unity in absolute value (we will 
conservatively assume that — 10 < Q < 10). As an illus- 
tration of the range of values allowed in some represen- 
tative models, Kaluza-Klcin-typc theories typically have 
1 < Q < 3, Einstein- Yang-Mills has Q = 1, and Randall- 
Sundrum type models have very small positive Qs (say 
Q ~ 0.01). These examples are discussed in more detail 
in |^|, [2^] — note that all of them have Q > 0. However, 
one can equally easily find models with Q < 
pie string theory dilaton-type models have Q 
while the BSBM-Brans-Dicke model has Q = 
0- 

In the rest of the paper we briefly sketch out anal ysis 
pipeline in Sec. [n] and present our results in Sec. HI. 
Finally we will discuss the implications of our results and 
present some conclusions. We emphasize that the key 
assumption of the present analysis is that the particle 
masses are kept constant, and there are no changes to 
the strong sector. A more general analysis, without this 
assumption, will be presented elsewhere. 
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while we constrain variations in the fine structure con- 
stant by considering the dimensionless parameter a/ao 
where ao is the fine structure constant today. 

The analysis method we adopt is based on the publicly 
available Markov Chain Monte Carlo package cosmomc 
PH with a convergence diagnostics done through the Gel- 
man and Rubin statistics. 

We sample the following ten-dimensional set of cosmo- 
logical parameters, adopting flat priors on each of them: 
the baryon and cold dark matter densities u>h and w c , 
the Hubble constant Ho, the scalar spectral index n s , 
the overall normalization of the spectrum A s at k = 0.05 
Mpc -1 , the optical depth to reionization, r and, finally, 
the variations in the fine structure constant a/cto and 
in the Newton's constant Xq- Furthermore, we consider 
purely adiabatic initial conditions and we impose spatial 
flatness. 

Our basic data set is the five-year WMAP data |32|. 
(temperature and polarization) with the routine for com- 
puting the likelihood supplied by the WMAP team. In 
addition to the WMAP data we also consider the fol- 
lowing CMB datasets: ACBAR (@), QUAD (gj|) 
and BICEP ( Pqj), as well as the older datasets from 
BOOMERanG (j|§) and CBI (§§]). For all these ex- 
periments we marginalize over a possible contamina- 
tion from Sunyaev-Zeldovich component, rescaling the 
WMAP template at the corresponding experimental fre- 
quencies. 

In what follows, we also combine the CMB data with 
the recent UNION catalog of supcrnovae type la lumi- 
nosity distances and with the improved constraint on the 
Hubble constant of h = 0.747 ± 0.036 at 68% c.L. from 
the recent analysis of (39). 

Constraints on Ag and a are also computed using stan- 
dard BBN theoretical predictions as provided by the nu- 
merical code described in [Q, which includes a full up- 
dating of all rates entering the nuclear chain based on 
the most recent experimental results on nuclear cross 
sections. The BBN predictions are compared with the 
experimental determinations of the 4 He mass fraction Y p 
and D/H abundance ratio, as discussed in |4l| 



0.250 ±0.003 



D/H = (2.87±g;l?) • 10" 



(3) 
(4) 



II. ANALYSIS METHOD 

We allow for a possible variation in the fine structure 
constant and in the Newton's constant during recombi- 
nation using the method described in [|l3] and |26| , mod- 
ifying the publicly available RECFAST (|§) routine in 
the CAMB (Q) CMB code. As in @ we consider vari- 
ations in the Newton's constant G by introducing a new 
dimensionless parameter Ag such that 



III. RESULTS 

In Table || we report the constraints on the a/ao and 
the Ag parameters obtained from the COSMOMC anal- 
ysis, using the the different combinations of the datasets 
described in the previous section, and in Figure |l| we 
show the 68% and 95% c.l. constraints on the a/ao vs 
Ac for the different datasets. 



3 



Experiment a/ao 68% c.l. Ag 68% c.l. 



All CMB 


0.999 


±0.017 


1.04 


±0.12 


All CMB+SN-Ia 


0.989 


±0.012 


1.04 


±0.11 


All CMB+HST 


1.003 


±0.008 


1.13 


±0.09 


ALL CMB+BBN 


0.985 


±0.009 


1.01 


±0.01 


Planck only 


1.000 


±0.015 


1.02 


±0.09 



TABLE I: Limits on a/ao and Ag from CMB data only (first 
row), from CMB+SN-Ia (second row), from CMB plus the 
HST prior on the Hubble constant, h = 0.748 ± 0.036 (third 
row), from CMB plus BBN (fourth row) and for simulated 
mock data for the Planck experiment. We report errors at 
68% confidence level. 
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FIG. 2: 68% and 95% c.l. constraints on the a/ao vs f° r 
different datasets. The contours regions come from CMB data 
(blue), CMB data and SN-Ia (red), and CMB±HST (green). 



Comparing these with the results of the recent studies 
for the two individual parameters [^l], |27j we see only a 
mild changes in the best fit and confidence intervals for 
a, while the changes are somewhat larger for G. 

More importantly, a degeneracy is clearly present be- 
tween Ag and the fine structure constant. The underly- 
ing reason is easy to understand. A change in a shifts 
the recombination epoch, affecting the angular diameter 
distance at recombination and the peaks position in the 
CMB anisotropy angular spectra. A similar effect can be 
obtained by changing the value of A<^ and the two pa- 
rameters are therefore degenerate. Both parameters are 
degenerate with the Hubble constant Ho as we can see 
from Figures || and [| 

From the above it follows that including the recent 
HST measurements of Hq has an important effect: it 
breaks the a-H n and Xq-Hq degeneracies and thereby 
provides a stronger bound on those parameters. This is 
clearly shown in the third row of Table 1 . 

1.1| 1 
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FIG. 1: 68% and 95% c.l. constraints on the a/ao vs Ag for 
different datasets. The contours regions come from CMB data 
(blue), CMB data and SN-Ia (red), and CMB±HST (green). 

As already mentioned in the introduction, from the 
perspective of mode-building it is interesting to consider 
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FIG. 3: 68% and 95% c.l. constraints on the Ag vs Ho for 
different datasets. The contours regions come from CMB data 
(blue), CMB data and SN-Ia (red), and CMB±HST (green). 

the behavior of the Q parameter, which in terms of our 
analysis pipeline is defined as 

--1 = Q(A 2 G -1). (5) 

Obviously, since the current data doesn't show any ev- 
idence for variations in G or a it is not possible to rule 
out any of the models discussed before. The CMB data, 
however, does show a clear correlation between a and Xq 
that could be fitted with a linear relation. Considering 
the models that are inside the 95% confidence level we 
found constraints on the best-fit Q parameter as reported 
in Table [i]. 

The data therefore prefers a value of Q whose value is 
significantly larger than zero, and again we see that the 
HST prior has a very strong effect. Taken at face value, 
this result would therefore strongly disfavor the models 
we mentioned in the introduction as having Q negative 
for Xq ^0.1 in absolute value. In other words, if there 
are any variations at or around the percent level, then 
the variations of the two parameters must have the same 
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Dataset Constraint on Q (at 68% c.l.) 



All CMB 


0.844 


< C 


} < 


0.888 


All CMB+SN-Ia 


0.868 


< c 


}< 


0.912 


All CMB+ HST 


0.744 


< c 


}< 


0.780 


Planck only 


0.872 


< c 


?< 


0.900 



TABLE II: Limits on the Q parameter from CMB data only 
(first row), from CMB+SN-Ia (second row), and from CMB 
plus the HST prior on the Hubble constant, h = 0.748 ±0.036 
(third row). We report errors at 68% confidence level. This 
analysis included only the models that are inside the 95% 
confidence level. 



sign. 

It is also interesting to notice that, due to the degen- 
eracy between the two parameters, a future detection for 
a variation in a could be on the contrary due to a varia- 
tion in G. It is therefore important to pursue a combined 
search for variations in the two constants since their ef- 
fect on the CMB anisotropy are very similar. In this 
respect we forecast the constraints achievable from the 
Planck satellite mission by using the specifications de- 
scribed in [^7], assuming a standard A-CDM model. 
These results are also presented in Tables I and II. As we 
can see, due to the a-G degeneracy the bounds obtained 
by Planck will be only marginally improved respect from 
current CMB data, as also shown in Figure 0. 




0.94 0.96 0.98, 1 1.02 1.04 



FIG. 4: 68% and 95% c.l. constraints on the Xg vs a/oto ex- 
pected from the Planck satellite alone (Blue) compared with 
the constraints obtained from current CMB data (Red). The 
degeneracy between these two parameters doesn't allow for 
significant improvements respect to current bounds. 



to variations of [i which may again have an imprint on 
the CMB, so a full analysis needs to be done before any 
stronger statements can be made. One analysis along 
these lines has been recently carried out j^], but the 
simplifying assumptions made by the authors imply that 
it applies only to a specific class of models, so a general 
analysis is still missing. 

Moreover, our analysis is made under the approxima- 
tion of a and G as constant in time and space. While the 
CMB constraints come mainly from recombination epoch 
and on scales larger than lOMpc, it is possible that our 
constraints could vary when a more accurate evolution in 
time and space up to the current epoch is considered. 

If we assume that a and G do not vary from BBN to 
recombination we can combine the CMB results with a 
BBN analysis. (This is a reasonable assumption, since in 
many scalar-field based models the field is frozen during 
the radiation era.) The results are plotted in Figure 5 
and also reported in Table 1. 

As we can see a sort of cosmic complementarity is 
present between the two datasets and much stronger 
bounds can be achieved. This can be understood as fol- 
lows. Differently than for the CMB, in the case of BBN 
in fact, variations of a and G are negatively correlated, 
since both Y p and Deuterium arc increasing functions of 
both parameters, see e.g. [Ell [h|. This implies that the 
likelihood contours for BBN and CMB are almost orthog- 
onal in the a — Xg plane, thus leading to a tighter bound, 
in particular on Xg- 
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FIG. 5: 68% and 95% c.l. constraints on the Xg vs q/qq ob- 
tained from current CMB data (Blue) compared with the con- 
straints obtained from current CMB data plus BBN bounds 
(Red). A cosmic complementarity between those two datasets 
helps in breaking the degeneracy and a much stronger bound 
is achievable on Xq- 



However, we again emphasize that these results were 
obtained on the assumption that a and G are the only 
rolling couplings, while everything else is kept fixed. 
While some of the above models should be taken as toy 
models and have little to say about this assumption, 
there are certainly others where one does expect par- 
ticle masses to vary. Among other effects, this will lead 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the effects of simul- 
taneous variations of the fine-structure constant a and 
Newton's gravitational constant G on the CMB. Our re- 
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suits indicate that the current data gives no clear indica- 
tion about the relative sign of the variations, but already 
prefers that any relative variations in a should be of the 
same sign of G for ~ 1% variations. We have also shown 
that, under the assumption that there's no rolling dur- 
ing the radiation era, much tighter constraints can be 
obtained by adding BBN data. 

While our analysis is more general (and more robust) 
than previous studies that have considered only the vari- 
ation of one of these couplings, it is by no means the final 
word on the subject. As a convenient simplification, we 
have assumed all particle masses to be constant. This 
assumption is often made (implicitly, if nothing else) in 
varying G studies since it gives an unambiguous meaning 
to its variations, but it is not a natural one. In particu- 
lar, one expects that in any sensible model where there 
are a variations, other quantities such as the proton-to- 
electron mass ratio [i will also vary. A more general anal- 
ysis, allowing for this possibility, is obviously much more 
difficult to carry out in full generality, but nevertheless 



still feasible. One such analysis will be presented in a 
follow-up paper. 

In any case, our results already show how fairly stan- 
dard astrophysical observablcs can place strong con- 
straints on high-energy physics models that would other- 
wise be difficult to test in laboratory or accelerator set- 
tings. With the significant gains in sensitivity expected 
for the next generation of ground and space experiments 
these constraints will become much stronger, and the 
early universe will become a key laboratory in which to 
probe fundamental physics. 



Acknowledgments 

The work of CM. is funded by a Ciencia2007 Re- 
search Contract, funded by FCT/MCTES (Portugal) and 
POPH/FSE (EC). 



[1] J. P. Uzan, Rev. Mod. Phys. 75 (2003) 403 [arXiv:hep- 
ph/0205340]; 

[2] E. Garcia-Berro, J. Isern and Y. A. Kubyshin, Astron. 

Astrophys. Rev. 14, 113 (2007). 
[3] M. T. Murphy, J. K. Webb and V. V. Flambaum, Mon. 

Not. Roy. Astron. Soc. 345, 609 (2003) [arXiv:astro- 

ph/0306483]. 

[4] E. Reinhold, R. Buning, U. Hollenstein, A. Ivanchik, 
P. Petitjean and W. Ubachs, Phys. Rev. Lett. 96, 151101 
(2006). 

[5] P. Molaro, S. A. Levshakov and M. G. Kozlov, Nucl. 
Phys. Proc. Suppl. 194, 287 (2009) [arXiv:0907.1192 
[astro-ph.CO]]. 

[6] R. Srianand, H. Chand, P. Petitjean and B. Aracil, Phys. 

Rev. Lett. 92, 121302 (2004) [arXiv:astro-ph/0402177]. 
[7] S. A. Levshakov et ai, Astronomy & Astrophysics 466, 

1077 (2007). 

[8] J. A. King, J. K. Webb, M. T. Murphy and 

R. F. Carswell, Phys. Rev. Lett. 101, 251304 (2008) 

[arXiv:0807.4366 [astro-ph]]. 
[9] R. I. Thompson et ai, Astrophys. J. 703, 1648 (2009) 

[arXiv:0907.4392 [astro-ph.CO]]. 
[10] T. Damour, F. Piazza and G. Veneziano, Phys. Rev. D 

66, 046007 (2002) [arXiv:hep-th/0205111]. 
[11] J. D. Barrow, J. Magueijo and H. B. Sandvik, Phys. Lett. 

B 541, 201 (2002) [arXiv:astro-ph/0204357]. 
[12] M. Nakashima, K. Ichikawa, R. Nagata and J. Yokoyama, 

arXiv:0910.0742 [astro-ph.CO]. 
[13] P. P. Avelino et ai, Phys. Rev. D 64 (2001) 103505 

[arXiv:astro-ph/0102144]; C. J. A. Martins, A. Mel- 

chiorri, R. Trotta, R. Bean, G. Rocha, P. P. Avelino 

and P. T. P. Viana, Phys. Rev. D 66 (2002) 023505 

[arXiv:astro-ph/0203149]. 
[14] S. Hannestad, Phys. Rev. D 60 (1999) 023515 

[arXiv:astro-ph/9810102]. 
[15] M. Kaplinghat, R. J. Scherrer and M. S. Turner, Phys. 

Rev. D 60 (1999) 023516 [arXiv:astro-ph/9810133]. 
[16] R. A. Battye, R. Crittenden and J. Weller, Phys. Rev. D 



63 (2001) 043505 [arXiv:astro-ph/0008265]. 

[17] C. J. A. Martins, A. Melchiorri, G. Rocha, R. Trotta, 
P. P. Avelino and P. T. P. Viana, Phys. Lett. B 585, 29 
(2004) [arXiv:astro-ph/0302295]; G. Rocha, R. Trotta, 
C. J. A. Martins, A. Melchiorri, P. P. Avelino, R. Bean 
and P. T. P. Viana, Mon. Not. Roy. Astron. Soc. 352, 20 
(2004) [arXiv:astro-ph/0309211]. 

[18] K. Ichikawa, T. Kanzaki and M. Kawasaki, Phys. Rev. D 
74 (2006) 023515 [arXiv:astro-ph/0602577]. 

[19] P. Stefanescu, New Astron. 12 (2007) 635 
[arXiv:0707.0190 [astro-ph]]. 

[20] M. Nakashima, R. Nagata and J. Yokoyama, Prog. 
Theor. Phys. 120 (2008) 1207 [arXiv:0810.1098 [astro- 
ph]]. 

[21] E. Menegoni, S. Galli, J. G. Bartlett, C. J. A. Mar- 
tins and A. Melchiorri, Phys. Rev. D 80, 087302 (2009) 
[arXiv:0909.3584 [astro-ph.CO]]. 

[22] P. R. Heyl, P. Chrzanowski, J. Res. Natl. Bur. Std. U.S. 
29, 1 (1942). 

[23] J. B. Fixler; G. T. Foster; J. M. McGuirk; M. A. Kase- 

vich, Science 315 (5808): 7477, 2007. 
[24] J. H. Gundlach, S. M. Merkowitz, Phys. Rev. Lett. 85, 

2869 (2000); T. J. Quinn, C. C. Speake, S. J. Richman, R. 

S. Davis, A. Picard, Phys. Rev. Lett. 87, 111101 (2001); 

W. Michaelis, H. Haars, R. Augustin, Metrologia 32, 267 

(1995) 

[25] O.G. Benvenuto, E. Garcia-Berro, and J. Isern, Phys. 
Rev. D 69 (2004) 082002. 

[26] O. Zahn and M. Zaldarriaga, Phys. Rev. D 67 (2003) 
063002 [arXiv:astro-ph/0212360]. 

[27] S. Galli, A. Melchiorri, G. F. Smoot and O. Zahn, 
Phys. Rev. D 80 (2009) 023508 [arXiv:0905.1808 [astro- 
ph.CO]]. 

[28] P. Loren-Aguilar, E. Garcia-Berro, J. Isern and 
Yu. A. Kubyshin, Class. Quant. Grav. 20, 3885 (2003) 
[arXiv:astro-ph/0309722]. 

[29] W. Y. Wong, A. Moss and D. Scott, arXiv:0711.1357 
[astro-ph] . 



6 



[30] A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 

538 (2000) 473 [arXiv:astro-ph/9911177]. 
[31] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002) 

(Available from http://cosmologist.info.) 
[32] G. Hinshaw et al. [WMAP Collaboration], 

arXiv:0803.0732 [astro-ph]. 
[33] E. Komatsu et al, arXiv:0803.0547 [astro-ph]. 
[34] C. L. Reichardt et al, Astrophys. J. 694 (2009) 1200 

[arXiv:0801.1491 [astro-ph]]. 
[35] M. L. Brown et al. [QUaD collaboration], 

arXiv:0906.1003 [astro-ph.CO]. 
[36] H. C. Chiang et al, arXiv:0906.1181 [astro-ph.CO]. 
[37] W. C. Jones et al, arXiv:astro-ph/0507494; F. Piacentini 

et al, arXiv:astro-ph/0507507; arXiv:astro-ph/0507514. 



[38] A. C. S. Readhead et al, Astrophys. J. 609, 498 (2004). 
[39] A. G. Riess et al, Astrophys. J. Suppl. 183 (2009) 109 

[arXiv:0905.0697 [astro-ph.CO]]. 
[40] O. Pisanti et al, Comput. Phys. Commun. 178 (2008) 

956 [arXiv:0705.0290 [astro-ph]]. 
[41] P.D. Serpico et al, 2004 JCAP 0412 010 [astro- 

ph/0408076]. 

[42] S. Galli, F. Iocco, G. Bertone and A. Melchiorri, 
Phys. Rev. D 80 (2009) 023505 [arXiv:0905.0003 [astro- 
ph.CO]]. 

[43] T. Dent, S. Stern and C. Wetterich, Phys. Rev. D 76 
(2007) 063513 [arXiv:0705.0696 [astro-ph]]. 



